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Abstract. Let (Z^' )tg[o,i] be a Hermite processes of order q and with 



Hurst parameter H £ {■^, 1). This process is H-self-similar, it has stationary 
increments and it exhibits long-range dependence. This class contains the 
fractional Brownian motion (for g = 1) and the Rosenblatt process (for q = 2). 
We study in this paper the variations of Z^i'^' by using multiple Wiener -Ito 
stochastic integrals and Malliavin calculus. We prove a reproduction property 
for this class of processes in the sense that the terms appearing in the chaotic 
decomposition of the their variations give birth to other Hermite processes of 
different orders and with different Hurst parameters. We apply our results to 

(^ , construct a consistent estimator for the self-similarity parameter from discrete 

-4-J ■ observations of a Hermite process. 



1. Introduction 



1.1. Background and motivation. The variations of a stochastic process play 

a crucial role in its probabilistic and statistical analysis. Best known is the qua- 

^O , dratic variation of a semimartingale, whose role is crucial in Ito's formula for 

seniiniartingales; this example also has a direct utility in practice, in estimating 
r — ^ , unknown parameters, such as volatility in financial models, in the so-called "his- 

r^ ' torical" context. For self-similar stochastic processes the study of their variations 

OO . constitutes a fundamental tool to construct good estimators for the self-similarity 

>^ ' parameter. These processes are well suited to model various phenomena where 

long memory is an important factor (internet traffic, hydrology, econometrics, 
among others). The most important modeling task is then to determine or esti- 
mate the self-similarity parameter, because it is also typically responsible for the 
C^ ' process's long memory and other regularity properties. Consequently, estimating 

this parameter represents an important research direction in theory and practice. 
Several approaches, such as wavelets, variations, maximum likelihood methods, 
have been proposed. We refer to the monograph [1] for a complete exposition. 
The family of Gaussian processes known as fractional Brownian motion (fBm) is 
particularly interesting, and most popular among self-similar processes, because 
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of fBni's stationary increments, its clear similarities and differences with standard 
Brownian motion, and the fact that its self-similarity parameter iJ, known as 
the Hurst parameter, is also clearly interpreted as the memory length parameter 
(the correlation of unit length increments n time units apart decays slowly at the 
speed n^^^^), and the regularity parameter (fBm is a-Holder continuous on any 
bounded time interval for any a < H). Soon after fBm's inception, the study of its 
variations began in the 1970's and 1980's; of interest to us in the present article are 
several such studies of variations which uncovered a generalization of fBm to non- 
Gaussian processes known as the Rosenblatt process and other Hermite processes: 
[2], [5], [7], [17] or [18]. We briefly recall some relevant basic facts. We consider 
(i?^)tg[o^i] a fractional Brownian motion with Hurst parameter H G (0,1). As 
such, B^ is the continuous centered Gaussian process with covariance function 

RH{t,s) =. E [Bf Bf ] ^ i(.2^ +t2^ - \t - s\^"), s,te [0, 1]. 

Equivalently, Bq ~ and E (^Bf — B/^) = \t — s\ . It is the only Gaussian 

process H which is self-similar with stationary increments. Consider = to < 
ii < . . . < <Ar = 1 a partition of the interval [0, 1] with t^ = -^ for « = 0, . . . , A^ 
and define the following q variations 
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where Hq with q > 2 represents the Hermite polynomial of degree q. Then it 
follows from the above references that for: 

• forO<-ff<l — ^ the limit in distribution of N^^v]^' is a centered 
Gaussian random variable, 

• for 7J = 1 — ^ the limit in distribution of {N log N)^^Vj^ is a centered 
Gaussian random variable, 

• for l — j-<H<l the limit oi N'^^^^-^^^^v)^ is a Hermite random variable 
of order q with self-similarity parameter q(H — 1) -I- 1. 

This latter random variable is non-Gaussian; it equals the value at time 1 of a 
Hermite process, which is a stochastic process in the qth Wiener chaos with the 
same covariance structure as fBm; as such, it has stationary increments and shares 
the same self-similarity, regularity, and long memory properties as fBm; see Def- 
inition 2.1. We also mention that very recently, various interesting results have 
been proven for weighted power variations of stochastic processes such as frac- 
tional Brownian motion (see [11]), fractional Brownian sheets (see [15]), iterated 
Brownian motion (see [12]) or the solution of the stochastic heat equation (see [16] 
or [3]). Because of a natural coupling, the last limit above also holds in L^(il) 
(see [13]). In the critical case H = \ — -^ the limit is still Gaussian but the nor- 
malization involves a logarithm. These results are widely applied to estimation 
problems; to avoid the barrier H = j that occurs in the case q — 2, one can use 
"higher-order filters" , which means that the increments of the fBm are replaced 
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by higher-order increments, such as discrete versions of higher-order derivatives, 
in order to obtain a Gaussian hniit for any H (see [9] , [8] , [4] ) . The appearance of 
Hermite random variables in the above non-central limit theorems begs the study 
of Hermite processes as such. Their practical aspects are striking: they provide 
a wide class of processes from which to model long memory, self-similarity, and 
Holder-regularity, allowing significant deviation from fBm and other Gaussian pro- 
cesses, without having to invoke non-linear stochastic differential equations based 
on fBm, and the notorious issues associated with them (see [14]). Just as in the 
case of fBm, the estimation of the Hermite process's parameter H is crucial for 
proper modeling; to our knowledge it has not been treated in the literature. We 
choose to tackle this issue by using variations methods, to find out how the above 
central and non-central limit theorems fit in a larger picture. 

1.2. Main results: summary and discussion. In this article, we show results 
that are interesting from a theoretical viewpoint, such as the reproduction prop- 
erties (the variations of Hermite processes give birth to other Hermite processes) ; 
and we provide an application to parameter estimation, in which care is taken to 
show that the estimators can be evaluated practically. Let Z'^'^'^^ be a Hermite 
process of order q with selfsimilarity parameter H € (^,1) as defined in Definition 
2.1. Define the centered quadratic variation statistic 
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(1.1) 



Also for H G (1/2,1), and q E N\ {0}, we define a constant which will recur 
throughout this article: 

{q\{H'{2H' -l))if'^ q 

We prove that, under suitable normalization, this sequence converges in L^(0) to 
a Rosenblatt random variable. 

Theorem 1.1. Let H e {1/2,1) andq eN\{0}. Let Z^'^^^^ be a Hermite process 
of order q and self- similarity parameter H (see Definition 2.1). Let c he an explicit 
positive constant (it is defined in Proposition 3.1). Then cN^'^^'^^'i'^Vn converges 



in L^ (Q.) to a standard Rosenblatt random variable with parameter H" := 

1, that is, to the value at time 1 of a Hermite process of order 2 and self- similarity 

parameter H" . 

The Rosenblatt random variable is a double integral with respect to the same 
Wiener process used to define the Hermite process; it is thus an element of the 
second Wiener chaos. Our result shows that fBm is the only Hermite process 
for which there exists a range of parameters allowing normal convergence of the 
quadratic variation, while for all other Hermite processes, convergence to a second 
chaos random variable is universal. Our proofs are based on chaos expansions 
into multiple Wiener integrals and Malliavin calculus. The main line of the proof 
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is as follows: since the variable Z^ is an element of the qth Wiener chaos, 

the product formula for multiple integrals implies that the statistics Vn can be 
decomposed into a sum of multiple integrals from the order 2 to the order 2q. The 
dominant term in this decomposition, which gives the final renormalization order 
j\j{2-2H)/q ^ is the term which is a double Wiener integral, and one proves it always 
converges to a Rosenblatt random variable; all other terms are of much lower 
orders, which is why the only remaining term, after renormalization, converges to 
a second chaos random variable. The difference with the fBm case comes from 
the limit of the term of order 2, which in that case is sometimes Gaussian and 
sometimes Rosenblatt-distributed, depending on the value of H. We also study 
the limits of the other terms in the decomposition of Vn, those of order higher 
than 2, and we obtain interesting facts: all these terms, except the term of highest 
order 2q, have limits which are Hermite random variables of various orders and self- 
similarity parameters. We call this the reproduction property for Hermite processes, 
because from one Hermite process of order q, one can reconstruct other Hermite 
processes of all lower orders. The exception to this rule is that the normalized 
term of highest order 2q converges to a Hermite r.v. of order 2q ii H > 3/4, but 
converges to Gaussian limit if i/ G (1/2, 3/4]. Summarizing, we have the following. 

Theorem 1.2. Let Z^'^'^> be again a Hermite process, as in the previous theorem. 
Let T2n be the term of order 2n in the Wiener chaos expansion of Vn : this is 
a multiple Wiener integral of order 2n, and we write Vn — X]n=i C2n72n where 
C2q-2k = kldY . 

• For every H € (1/2,1) and for every k — l,...,g — 1 we have con- 
vergence of the expression {zk,H)~ TV^^"^^ •"^T'2fe in L^ (fi) to Z'-^-'^\ a 
Hermite random variable of order r ~ 2k with self- similarity parameter 
K = 2k{H' — 1) + 1, where Zk.H is a constant. 

• For every H G (1/2, 3/4), with k — q, we have convergence of x^ ^ „\NT2q 
to a standard normal distribution, with Xi^H,q o, positive constant depending 
on H and q. 

• For every H G (3/4, 1), with k = q, we have convergence of X2 i^ „ N^^^^ 
T2q in L^(il) to ^2g,2_f/-i^ ^ Hermite r.v. of order 2q with parameter 
2H — 1; with X2.H,q o. positive constant depending on H and q. 

• For H = 3/4, with k = q, we have convergence of y/N/ log Nx^ ^ T2q to 
a standard normal distribution, with x^h o, positive constant depending 
on H and q. 

Some of the aspects of this theorem had been discovered in the case oi q = 2 
(Rosenblatt process) in [20]. In that paper, the existence of a higher-chaos-order 
term with normal convergence had been discovered for the Rosenblatt process with 
H < 3/4, while the case oi H > 3/4 had not been studied. The entire spectrum 
of convergences in Theorem 1.2 was not apparent in [20], however, because it 
was unclear whether the term r2's convergence to a Rosenblatt r.v. was due 
to the fact that we were dealing with input coming from a Rosenblatt process, 
or whether it was a more general function of the structure of the second Wiener 
chaos; here we see that the second alternative is true. The paper [20] also exhibited 
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a remarkable structure of the Rosenblatt data when H < 3/4. In that case, as 
we see in Theorem 1.2, there are only two terms in the expansion of Vjv, T2 and 
T4; moreover, and this is the remarkable feature, the proper normalization of 
the term T2 converges to none other than the Rosenblatt r.v. at time 1. Since 
this value is part of the observed data, one can subtract it to take advantage of 
the Gaussian limit of the renormalized T4, including an application to parameter 
estimation in [20]. In Theorem 1.2 above, ii q > 2, even li H < 3/4, by which 
a Gaussian limit can be constructed from the renormalized T2q, we still have 
at least two other terms T2,T4, • • • T2g_2, and all but at most one of these will 
converge in L^ (fl) to Hermite processes with different orders, all different from 
q, which implies that they are not directly observed. This means our Theorem 
1.2 proves that the operation performed with Rosenblatt data, subtracting an 
observed quantity to isolate T2q and its Gaussian asymptotics, is not possible with 
any higher-order Hermite processes. The last aspect of this paper applies Theorem 
1.1 to estimating the parameter H. Let Sn be the empirical mean of the individual 
squared increments 

and let 

HN = ilogSN)/i2\ogN). 

We show that Hn is a strongly consistent estimator of H, and we show asymp- 
totic Rosenblatt distribution for N'^C^-"^/'^ (h - h") log TV. The fact that this 
estimator fails to be asymptotically normal is not a problem in itself. What is 
more problematic is the fact that if one tries to check ones assumptions on the 
data by comparing the asymptotics of H with a Rosenblatt distribution, one has 
to know something about the normalization constant N^^^'^yi, Here, we prove 
in addition that one may replace H in this formula by Hj^, so that the asymptotic 
properties of Hn can actually be checked. 

Theorem 1.3. The estimator H^ is strongly consistent, i.e. limjv-s-oo -ffw = H 
almost surely. Moreover there exists a standard Rosenblatt random variable R with 
self-similarity parameter 1 + 2{H — l)/q such that 



lim E 



Here a {H') = {I + {H - l)/q) (1 + 2{H - l)/q) . 



= 0, 



Replacing the constant ci,h by its value in terms of H^q instead of H also seems 
to lead to the above convergence. However, this article does not present any nu- 
merical results illustrating model validation based on the above theorem; moreover 
such applications would be much more sensitive to the convergence speed than to 
the actual constants; therefore we omit the proof of this further improvement on 

Cl, H. 
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2. Preliminaries 

2.1. Multiplication in Wiener Chaos. Let {Wt)telo.i] be a classical Wiener 
process on a standard Wiener space {il,T,P). If / G i^([0,l]'") with m > 1 
integer, we introduce the multiple Wiener-Ito integral of / with respect to W. 
We refer to [14] for a detailed exposition of the construction and the properties 
of multiple Wiener-Ito integrals. Let / G S, i.e. / is an elementary function, 
meaning that 

./ ^ / ^ Cil,...im l^i. X...X/li^ 

ii .....ijn 

where ii, . . . ,im describes a finite set and the coefficients satisfy c^j^ ^^ — if two 
indices ik and ii are equal and the sets Ai G B{[Q, 1]) arc disjoints. For such a step 
function / we define 

Im{f)= J2 C^„...^^W{A^,)...WiA^J 

where we put iy([a, 6]) — Wb — Wa- It can be seen that the appfication /^ 
constructed above from S to L^(i7) satisfies 

E[/„(,/)/™(.9)] =n!(/,5)L^([o,i]") ^ rn = n (2.1) 

and 

E [/„(/)/„ (.9)] =0 if m^n. 

It also holds that /„(/) = /„ (/) where / denotes the symmetrization of / de- 
fined by /(xi, . . . ,x„) = ^ J2ae„ fi^cr{i),- ■ ■,x^^n))- Since the set S is dense in 
L^([0, 1]") the mapping /„ can be extended to a linear continuous operator from 
from i^([0, 1]") to L^(ri) and the above properties hold true for this extension. 
Note also that /„(/) with / symmetric can be viewed as an iterated stochastic 
integral 

In{.f) = nl f r ... f' f{h,...,t„)dWt,...dWt„; 
Jo Jo Jo 

here the integrals are of Ito type; this formula is easy to show for elementary /'s, 

and follows for general symmetric function / G ^^([0, 1]") by a density argument. 

We recall the product for two multiple integrals (see [14]): if / G i^([0, 1]") and 

g G i^([0, 1]™) are symmetric then it holds that 

rnAn / \ / \ 
In{f)lm{9) = E ''(7) \,)lm+n-2l{f®l9) (2.2) 



(=0 



where the contraction f(E)ig belongs to i^([0, 1]™+" 2(-j for / = 0, 1, . . . , mAn and 
it is given by 

(/ "^l 9){S1, . . . , Sn-hh, . . . , tm-i) = 

/(si, . . . , s„_;. Ml, . . . , ui)g{ti,. . . , tm-hUi, . . . , ui)dui ...dui. (2.3) 

[0,1]' 
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2.2. The Hermite process. Recall that the fractional Brownian motion process 
{Bl^)te[o,i] with Hurst parameter H E (0, 1) can be written as 

5f ^ f K"{t,s)dWs, t e [0, 1] 
Jo 

where {Wt,t <E [0,T]) is a standard Wiener process, K^ {t,s) = c^s^^^ J {u — 

s)^^^u^^^du ii t > s (and it is zero otherwise), ch — ( b(2--2H^--^) I ^^"^ 

/?(•,•) is the Beta function. For t > s, the kernel's derivative is ^M^{t,s) = 

Ch (f)^ {t — s)^~2. Fortunately we will not need to use these expressions 
explicitly, since they will be involved below only in integrals whose expressions are 
known. 

We will denote by {Z^' )t^m^i] the Hermite process with self-similarity param- 
eter H € (1/2, 1). Here g > 1 is an integer. The Hermite process can be defined 
in two ways: as a multiple integral with respect to the standard Wiener process 
{Wt)te[a.i\\ 01' 8-3 a multiple integral with respect to a fractional Brownian motion 
with suitable Hurst parameter. We adopt the first approach throughout the pa- 
per. We refer to [13] of [19] for the following integral representation of Hermite 
processes. 

Definition 2.1. The Hermite process {Z^' )tG[o,i] of order q > 1 and with 
self-similarity parameter H £ (^ , 1) is given by 

Zl"'"^ ^d{H,q) f dWy,...dWy^ ( [ diK""' {u,y^) . . .diK"' {u,yq)du\ 

J[Q,t]i \JyiV...Vv, J 

(2.4) 
for t E [0, 1], where K^ is the usual kernel of the fractional Brownian motion and 

TT 1 

H' = 1 + -^=^ {2H' - 2)q = 2H -2. (2.5) 

Of fundamental importance is the fact that the covariance of Z'^'^^ is identical 
to that of fBm, namely 

E [zi^'^^Z^,^'"^] = -it'" + s2« - |t - 5|2ff). 

The constant d{H,q), given in (1.2) on page 3, is chosen to avoid any additional 
multiplicative constants. We stress that Z'"^'^' is far from Gaussian for g > 1, 
since it is formed of multiple Wiener integrals of order q. 

The basic properties of the Hermite process are listed below: i) the Hermite 
process Z^'^'"^ is iJ-selfsimilar and it has stationary increments; ii) the mean square 
of the increment is given by 

2" 



E 



Zi^,H) _ ^igM) 



s 



\2H. 



(2.6) 



as a consequence, it follows from Kolmogorov's continuity criterion that Z^'^'"^ has 
Holder-continuous paths of any order S < H; iii) it exhibits long-range dependence 
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in the sense that J2n>i^ ^i i^n+i ^ ^^' ) ^ °°- I^ f^*^t, the summand 
in this series is of order r?^"^ . This property is identical to that of fBm since the 
processes share the same covariance structure, and the property is weh-known for 
fBm with H > 1/2; iv) li q = 1 then Z^^'^' is a standard Brownian motion with 
Hurst parameter H while for q > 2 this stochastic process is not Gaussian. In the 
case q ^ 2 this stochastic process is known as the Rosenblatt process. 



3. Variations of the Hermite process 



Since E (z'f+"^ - z'-f"A ^ N-^" and by (2.6), the centered quadratic vari- 
ation statistic Vat given in the introduction can be written as 



N-l 



Vn = N^"-' J2 



^(qM) _ ^{q,H) 



N 



-2H 



Let li := \j^, ^^]- In preparation for calculating the variance of Vn we will find 
an explicit expansion of Vn in Wiener chaos. We have Z^+i — Z_f ' — Iq (/i,Ar) 
where we denoted by fi^N^yi, ■ ■ ■ , Vq) the expression 



-H' I 



^H' I 



l[o,i+i] (yi V . . . V yq)d{H, q) / diK"" {u, yi) . . . diK"" {u, yq)du 

' " JyiV...Vyq 



-H' I 



-H' I 



-l[o.^] (yi V . . . V yq)d{H, q) / diK"" {u, yi) . . . ^iX^ (u, yq)du. 



Using the product formula for multiple integrals (2.2), we obtain 

1 / \ 2 

Iq{ft.N)Iq{fi,N) =X!^'(fc) ^2g-2( {ft,N (^l Ji.n) 

where the f (E)ig denotes the Z-contraction of the functions / and g given by (2.3). 
Let us compute these contractions; for I = q we have 



(/i,Ar (Xlq fi,N) = q^-{fi,N , Ii.n) L^{[id,l])S>i = E 



\ N 



z 



iq.H) 



N' 



-2H 



Throughout the paper the notation diK{t,s) will be used for diK^ (t,s). For 
Z = we have 



{ft.N ^O fi,N)iyi ,...yq,Zl,...,Zq) = {fi^N «) fi,N)iyi , . . . y,, Zl , . . . , Zq) 

= fi,N{yi,---,yq)fi.N{zi,...,Zg) 
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while for 1 < k < q — 1 

{fi,N<^k fi.N){yi,---yq-k,Zi,...,Zq^k) = d{H,q)'^ / dai...dak 

J[0.1]'' 

li+i,g-fcli+i,fe / dudiK{u, yi)... diK{u, yg-k)diK{u, ai) . . . diK{u, ak) 

- lf;^_fel"fc / dudiK{u, yi)... diK{u, yg^k)diK{u, ai) . . . diK{u, au) 

i,k / 

^+i,q-k'^t+i,k / dvdiK{v, zi)... diK{v, Zq_k)diK{v, ai) . . . diK{v, ak) 

\ ''^Ui.k 

- ^Tq-k'^Tk I " dvdiK{v, zi) . . . diKiy, Zg.k}diK{v, ai) . . . diK{v, at) ] 



'''- ^ IV 



where 1^\ denotes the indicator function Itq _L^k(xj) with x being y, 2;, or a, and 
I^j, denotes the interval [xi V • • • V Xq-k V ai • • • V ak',i/N], with x being y or z. 
By interchanging the order of the integration we get 

{.fi,N ®k fi,N){yi, ■ ■ ■ yq-k, zi,..., Zg-k) 

= c?(-ff, 9) S l[o,i+l]2,-2fc(yi,2;i) / dudiK{u,yi) ...diK{u,yq^k) 

/ /^uAv \ ^ 

dvdiK{v, zi) . . .diK{v, Zg^k) ( / diK{u,a)diK{v,a)da\ 
ziv...zg_fc \Jo J 

-l[o,i+l]9-fc(yi)l[o,j^]'j-'=(^i) / dudiK{u,yi) . . .diK{u,yq^k) 

r-k / /-"^^ n'' 

/ dvdiK{v, zi) . . .diK(v, Zq^k) ( / diK{u,a)diK{v,a)da\ 

JziV...z,_fc \Jo ) 

-l[o,^]9-fc(yj)l[o,i+i]5-*(2i) / dudiK{u,yi) . . .diK{u,yg_k) 

/ dvdiK{v, zi) . . .diK{v,Zg^k) { / diK{u,a)diK{v,a)da\ 

JziV...z,_fc \Jo J 

+l[o,^]9-(= (y*)l[o,^]9-»= (zi) / dudiK{u, yi)... diK{u, y,_fe) 

•^ J/lV...J/<,_fc 

(ii;9iii'(v, zi) . . . 9i-ftr(u, Zg_fe) / diK{u,a)diK{v^a)da 
and since 

/^uAv 

.\2H'-2 



pU/\V 

I diK{u,a)diK{v,a)da — a{H')\u — 
Jo 



V 
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with a{H') ^ H'{2H' - 1), wc obtain 

{fi,N ®k fi,N) (yi, • • • yq-k,zi, ■■■, Zq-k) ^ d{H, qfa{H'f 



X < l[o_i+i],-fc (2/i)l[o,^]'!-'= i^i) I dudiK{u, j/i) . . . diK{u, Vq-k) 



dvdiK{v, zi) . . . diK{v, z,^k)\u - i;|(2ff'-2)fe 



Zl\/ ...Zg^k 



Jyi'^-'-Vq-k 



X / dvdiK{v^zi) . . .diK{v^Zq^k)\u — v\ 

'ziV...z„_fc 



{2H'-2)k 



-l\o^^],-k{yi)l 



[0,1+1] 9- 4^0 / dudiK{u,yi) . . .diK{u,yg-k) 

Jj/lV...j/,_fc 



X / dvdiK{v,zi)...diK{v,Zq^k)\u-v\'^^"'-^^'' 

I ZlV...Zg-k 

fir 

+ l[o,^],-'.(y*)l[o,^]9-fc(20 / dudiK{u,yi) . . .diK{u,yq^k) 



_„,U2H'-2)k 



X I dvdiK{v,zi) . . .diK{v,Zq^k)\u — v\ 

'ZlV...Zg_fc 



As a consequence, we can write 

Vat = T2q + C2q-2T2q-2 + • • • + C^T^ + C2T2 

where 



C2q-2k 



k\ 



(3.1) 

(3.2) 
(3.3) 



are the combinatorial constants from the product formula for < fc < g — 1, and 

T2q-2k ■■= N^"-'l2q-2k ( Y. A^ ^k /^^ I ' (3-4) 



. i=0 



where the integrands in the last formula above are given explicitly in (3.1). This 
Wiener chaos decomposition of Vn allows us to find V^'s precise order of magni- 
tude via its variance's asymptotics. 



Proposition 3.1. With 



Cl.H 



Ad{H,q)'^{H'{2H' -l)f<i 



{AH' - 3)(4H' - 2)[{2H' - 2){q - 1) + 1]^[{H' - l)iq - 1) + 1]2 ' 
it holds that 



(3.5) 



lim E 



-1 jsr(2-2H')2 -2y2 



= 1. 
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Proof. Wc only need to estimate the L^ norm of each term appearing in the 
chaos decomposition (3.2) of Vn , since these terms are orthogonal in L^ . We can 
write, for < fc < 9 — 1, 



^[Tl-2k] - N^"-\2q-2k)\ 



^N~l 



X! /»,JV ®k fi 



N 



i=0 
AT-l 



L2([0,l]29-2fc) 



N {2q—2k)\ y ^ {fi,N'Sikfi,N, fj,N'Sikfj,N)L'^(\oj]-^i-2k) 

ij=0 



where (gY — g and fi^N®kIi,N denotes the symmetrization of the function fi^N®k 
fi,N- We will consider first the term T2 obtained for k = q — 1. In this case, the 
kernel X^j^q f^-^ '^q-i fi.N is symmetric and we can avoid its symmetrization. 
Therefore 



N-l 



E [T^] - 2! N^"-^ J2 -f^^^ ®9-i /'^A^IIi^([o,i]^) 



i=0 

N-1 



2^- N 2_^ {fi,N (Xig-l fi,N,fj,N <^q-l /j,Af)L2([oa]2). 



We compute now the scalar product in the above expression. By using Fubini 
theorem, we end up with the following easier expression 



{ft^N -Siq-i fi.Njj.N ®q-i /j-.Ar)L2([o.il2) = a{H'f'id{H,qf 

Jii Jii Jij J ij 

X\U- ^,|(2ff'-2)(«-l) 1^,' _ ^'\(2H'-2)(q-l) |,^ _ ^'|2H'-2|^ _ v']''"' '^ dv' du' dvdu 



Using the change of variables y = {u— jf)N and similarly for the other variables, 
we now obtain 



E [r|] = 2d{H,qt{H'{2H' - l))29^4H-2^-4^-(2H'-2)2g 
A^-l ,.1 ,.1 ,.1 ,.1 

xE 



ij=0 



JO JO Ja 



dydzdy'dz'ly - z|(2«'-2)(«-i) |y' _ ^'|(2ff'-2)(,-i) 



:\y-y' + ^-J\^'"'-'^\z-z' + ^-J\^'"'-'\ 



This can be viewed as the sum of a diagonal part {i — j) and a non-diagonal part 
(i ^ j), where the non-diagonal part is dominant, as the reader will readily check. 
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Therefore, the behavior of E [r|] will be given by 
E [Tf ] 

:= 2\d{H,q)%H'{2H' -I)fm-^2J2 [ [ [ I dydzdy'dz' 

.^j Jo Jo Jo Jo 

N-2N-i „i „i „i „i 

= 2\diH,qf{H\2H'-l)fm-'2j2J2 / / dydzdy'dz' 

i=0 e=2 " '^ 

x{\y - z\\y' - z'|)(2^'-2)(«-i)(|y -y' + £-l\\z-z' + £- l|)(2^'-2) 

= 2\d{H,qf{H'{2H'-l)f'iN-^2'Y^{N-l+\) / / / dydzdy'dz' 

x(|y-z||y'-z'|)(2^'-2)(^-i)(|y-y' + ^-l||z-z' + £-l|)(2^'-2). 
As in [20] note that 

1 ^ 

w 



7V2 



£=2 



1 ^ 



e-l y-y' ^- l|2H'-2,^-2' ^-1 



2H'-2 



1=2 



N " N 



N 



N N 



Using a Riemann sum approximation argument we conclude that 



E 



T. 



'2 



^d{H,q)\H'{2H' -l)f<i X Ar2(2H'-2) 



(4ff' - 3)(4i/' - 2)[{{2H' -2){q-l) + l)Y[{H' - l){q - 1) + If ' 
Therefore, it follows that 



E 



-1 j^2(2-2H')rj.2 



>N^oo 1, 



(3.6) 



with ciH as in (3.5). 

Let us study now the term T/^, . . . , T2q given by (3.4). Here the function 

X]i=o /«.^ '^'^ /^^ ^'^ ^'^ longer symmetric but we will show that the behavior of 
its L^ norm is dominated by E [Tj^] . Since for any square integrable function g 
one has \\g\\L^ < \\g\\L^, we have for k = 0, . . . ,q — 2 



1 



N-l 



{2q-2k) 



E \Ti 



-9i-V^L^2q-2fc 



= A^' 



AH-2\ 



7 , /i,JV®fc/i,Ar|li2([o,l]2<!-2fc) 



i=0 



N-1 



< ^^""-'ll^ AjV^fc/, 



i,A^llL2([0,l]25-2fc) 



1=0 
JV-1 



N 2_^ {.fi^N (Xlfe /i,Ar, /j,JV (^Jfe /j-^Ar)i2([o,l]29-2fc) 



(3.7) 
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and proceeding as above, with eH,q,k '■= {'2q — 2k)l{H'{2H' — l))^'^d{H, g)^ we can 
write 



N-l 



E [T,V2.] < eH,,.,N*"'' E 



i,j=0 



dyidzi 



h -Ih 



/, J I 



dy'idz'i 



and using a change of variables as before, 
E [^"25-2*:] 



< 



N-l „ 

eH,,.kN^"-^N-^N-^^"'-^^^^ V / dydzdy'dzWy - z\\y' - z'\) 



I _ l\\{2H'-2)k 



y<\y-y' + l- j|(2ff'-2)(9-fe) \^-^' + ^- ^■|(2H'-2)(,-fe) 



(2g - 2k)\d{H, qY iv(2ff'-2)(29-2fc2 ^ , ^ 



aiH')-^i 



N^ 



£-1 

N 



x{\y^z\\y'^z'\r"'-'^' 



y-y' , ^-1 



[04]" 
^-1 



iV 



iV 



dydzdy' dz' 

(2H'-2){q-k) 

(3.8) 



TV TV 

Since off a diagonal term (again of lower order), the terms ■^^- are dominated by 
-^ for large l,N it follows that, for 1 < fc < g — 1 



E 



-q-kM 



iy{2-2H'){2q-2k)rp2 



2q-2k 



0(1) 



(3.9) 



when iV — > 00, with 



Cq-fe,H -2(/ (l-x)x(2^'-2)(29-2fc)(|^)a(^')-2(2^_2fc)!d(i7,(7)2a(H')29. (3.10) 
Jo 

It is obvious that the dominant term in the decomposition of Vn is the term in 
the chaos of order 2. [The case A: = is in the same situation for i? > | and for 
H e (i, |) the term T2q obtained for fc = has to be renormalized by A^, (as in 
proofs in the remainder of the paper); in any case it is dominated by the term T2]. 
More specifically we have for any k < q ~ 2, 



E 



A72(2-2H')rp2 

^^ ^2q-2k 



O f ]\j-^{^-^H')^(i-'''-^)\ 



(3.11) 



Combining this with the orthogonality of chaos integrals, we immediately get that, 
up to terms that tend to 0, N"^^"^^ Vn and 7V2~2ff j^^ have the same norm in 
L^ (f2). This finishes the proof of the proposition. ■ 



Summarizing the spirit of the above proof, to understand the behavior of the 
renormalized sequence Vn it suffices to study the limit of the term 



h iN^"-'N^^-^"'^Y.f^-.^®^~^f^A 



(3.12) 
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with 

{.h,N <^q-i fi,N){y, z) (3.13) 

(l[o,^](yVz) / / dvdudiK{u,y)diK{v,z)\u-v\^'^"'-'^^^'i-^^ 
+l[o,^l(y)l/.Wy J " dvdudiK{u,y)diK{v,z)\u-v\^^"'-^^('^-^'> 

i + l 

+l/,(y)l[o,^]W / " / dvdudiK{u,y)d,K{v,z)\u-v\(^"'-^^'^'i-'') 

Jy J li 

+l7,(y)l/,(z) / " / " dvdudiK(u,y)diK(v,z)\u^v\'^^"'-^^('^-^'>\. 

We will see in the proof of the next theorem that, of the contribution of the four 
terms on the right-hand side of (3.13), only the first one does not tend to in 
L^ (n). Hence the following notation will be useful: f^ will denote the integrand 
of the contribution to (3.12) corresponding to that first term, and r2 will be the 
remainder of the integrand in (3.12). In other words, 

7V-1 

/f + r2 = ^2ff-i^(2-2H') J2 /,,^ 0,_i /,,^ (3.14) 

and 



/f(y,z) : =N^"-^N(^-^"">d{H,qfa{H'Y~^ 
X]l[o,T^](yV2) /" f dvdudiK{u,y)diK{v,z)\u-v\^'^"'-'^^^'i-^\ (3.15) 



Af-l 

i[o,^](yVz) 

i=0 •'^i '' ^i 

Theorem 3.2. The sequence given by (3.12) converges in 1? (f2) as iV —> oo to the 
constant c-^^ times a standard Rosenblatt random variable Z\ ' with self- 

similarity parameter 2H' ~ 1 and H' is given by (2.5). Consequently, we also have 
that c^ Jj ]^{'^-'^H )^- y^ converges in L^ (Jl) as N ^^ co to the same Rosenblatt 
random variable. 

Proof: The first statement of the theorem is that N^^"^^ T2 converges to 

1/2 ^(2,2ff'-l) 

in L'^{il). From (3.12) it follows that T2 is a second-chaos random variable, with 
kernel 

N'""-' J2 ^f^-« ®«-i ^^■^) 

(see expression in ( 3.13)), so we only need to prove this kernel converges in 
L^ ([0,1]^). The first observation is that r2{y,z) defined in (3.14) converges to 
zero in L^([0, 1]^) as A^ -^ 00. The crucial fact is that the intervals li which are 
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disjoints, appear in the expression of this term and this implies that the non- 
diagonal terms vanish when we take the square norm of the sum; in fact it can 
easily be seen that the norm in L^ of r2 corresponds to the diagonal part in 
the evaluation in ET2 which is clearly dominated by the non-diagonal part, so 
this result comes as no surprise. The proof follows the lines of [20]. This shows 
j^(2-2H )rp^ jg |.]_^g g^j^ Qj: j^ {fi^) &n.d a term which goes to in L^ (Q). Our next 
step is thus simply to calculate the limit in L^ (57), if any, of /2(/2^) where f^ 
has been defined in (3.15). By the isometry property (2.1), limits of second-chaos 
r.v.'s in L^ (O) are equivalent to limits of their symmetric kernels in L^ ([0, 1]^). 
Note that f^ is symmetric. Therefore, it is sufficient to prove that f^ converges 
to the kernel of the Rosenblatt process at time 1. We have by definition 



f^{y,z) = {H'{2H' -l)f-^-'-^d{H,qfN^"-^N''-^''' 



X 



Thus for every y, z, 



^-^ f f 

^H \u-v\^^"-^^^''-^'>diK" {u,y)diK" {v,z)dudv 



N-l 

■ , _^|(2H'-2)(,-l) 



d{H,qf{H'{2H'-l))<^'^-'^N^"-'N^-^"' J2 

1=0 

X (diK"' (u, y)diK"' (v, z) - diK"' {i/N, z)diK"' {i/N, z) \ dudv 



+diH,qfiH'i2H'-l)Y''-^'>N^"-^N^-^"'J2l I l«-«l 

X diK"' {i/N, y)diK"' {i/N, z)dudv 
AUy,z)+A^{y,z). 



(2H'-2)(q-l) 



As in [20], one can show that E 



\aN\ 



as A^ ^ 00. Regarding the 



|2 

Il2([04]2)J 

second term A2 {y, z), the summand is zero if i/N < y V z, therefore we get that 
/^ is equivalent to 



N^"-^N^-^"'d{H,qf{H'{2H' - l))'"-!) Y. 

1=0 

xdiK"' {i/N,y)diK"' {i/N,z)dudv 



N-l 

(2H'-2)(q-l) 



\U — V\ 
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t2-2H' 



{H'{2H' - l))(«-i)d(i7,q)2iV2ff-i7v2 



N-1 



xJ2diK'''ii/N,y)d,K"'{i/N,z)ly^,^,/NJ j \u - v\'^^"'-^^^^-'Uudv 



4=0 



{H'{2H' - 1))(''-I)[((2i7' - 2)(q - 1) + mn' - l){q - 1) + 1)]-1 



d{H, qf 



J2 diK"' ii/N,y)diK"' ii/N,z)ly^,^,/N 
{H'{2H' - l))(«-i) 



d{2H' - 2, 2) {{2H' - 2){q - 1) + 1){{H' - l){q - 1) + 1) 



N-l 



id{2H' -2,2)N-^ ^ diK"' {i/N,y)diK"' {i/N,z)l 



yVz<i/N- 



i=0 



r-iY-JV-la T^H'u, 



The sequence d{2H' - 2,2)N-^J2"=a d^^ {i/N,y)diK" {i/N,z)ly^^<,fN is a 
Riemann sum that converges pointwise on [0,1]^ to the kernel of the Rosenblatt 
process Z^^ ~^-^ at time 1. To obtain the convergence in 1? ([0, 1]^) we will apply 
the dominated convergence theorem. Indeed, 



nl ,.1 -, ^-1 



Jo Jo ^^ ,_n 



/o Jo 

1 JV-1 



yVz<i/N 



iV2 



i,j=0 
N-1 



i=0 

1 



dydz 



'H' I ■ , 



d^K'^ {i/N,y)diK" (j7^,2/)ly<(4A,)/Arrfy 



^ 4E|E[A^V^^AZ,7^^]|^ 



i,j=0 



where AZ^^j.^ is the difference Z (-^) — Z (■^) for a Rosenblatt process Z. We now 
show that the above sum is always ^ iV^, which proves that the last expression, 
with the N~'^ factor, is bounded. In fact for Hi = 2H' — 1 



N-l N-1 

i,j— "ijj— 



i-j + l 



N 



2 Hi 



i — j — 1 



N 



2Hi 



i- J 



N 



2Hi 



N-l 



2Hi 



= -^J2 \\^-J + M'"'+\^-J-M'"'-'2\^-J\ 

i,j=0 
ivr-4Hi N-l 

< 27V J2 W + if"' + \i-if"' -'2\e\'^"' 

i=-N+l 



The function 5 (£) = \£+l\^"^ + \£-l\^"^ -2\£\^"^ behaves like i/i(2i7i-l)|£|2ffi-2 
for large £. We need to separate the cases of convergence and divergence of the 
series X^^oo l.9(^)P- I* i^ divergent as soon as Hi > 3/4, or equivalently H' > 7/8, 
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in which case we get for some constant c not dependent on N, 

N-l 

The series X^^^oo l5(^)P i^ convergent if H' < 7/8, in which case we get 

JV-l 

J2 |E [AZ,/nAZ^/n] r < cN-^"^+\ 

ij=0 

For this to be <C N'^, we simply need —AHi + 1 < 2, i.e. H' > 5/8. However, since 
q > 2 and H > 1/2 we always have H' > 3/4. Therefore in all cases, the sequence 
A^(j/, z) is bounded in L^ ([0, 1]^) and in this way we obtain the L"^ convergence 
to the kernel of a Rosenblatt process of order 1. The first statement of the theorem 
is proved. In order to show that c^ ^ ]\[i'^-'^H )q-^Yj^ converges in LF'iVl) to the 
same Rosenblatt random variable as the normalized version of the quantity in 
(3.12), it is sufficient to show that, after normalization by A^^^^h ^ each of the 
remaining terms of in the chaos expansion (3.2) of Vjv, converge to zero in L^(0), 
i.e. that A^(2-2_ff )T2q-2k converge to zero in L^(fi), for aU 1 < fc < q - 1. From 
(3.11) we have 



E 



^2(2-2H')rJ.2^_^^^ ^ Q |^^-2(2-2ff')2(9-fc-l)^ 



which is all that is needed, concluding the proof of the theorem. 



4. Reproduction property for the Hermite process 

We now study the limits of the other terms in the chaos expansion (3.2) of 
Viv. We will consider first the convergence of the term of greatest order T^q in 
this expansion. The behavior of T2q is interesting because it differs from the 
behavior of the all other terms: its suitable normalization possesses a Gaussian 
limit if i/ G (1/2,3/4]. Therefore it inherits, in some sense, the properties of the 
quadratic variations for the fractional Brownian motion (results in [20]). We have 



T2q = N^"-'l2q J2 -f^^^ ® -^^^ 






and E [Ti^] = N^"-^{2q)\Y.li=o{kN®kN, fj,N®fj,N)L-[o,iY^- We will use the 
following combinatorial formula: if f,g are two symmetric functions in L^([0, 1]"^), 



(2g)! (/«)/, g(8)5r)i2([o,i]2, 

9-1 



9-1 / \ 2 



k=l 
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to obtain 



N' 



4H-2 



N-1 



JV-1 g-l / \ 2 
-N ' (ql) 2^ 2^ ( , I {fi,N (i^k f],N,fjM (^k fi,N)L2{[0,l]2i-^'')- 



First note that E 



ij=a fc=i 



(qM) r,iq.H) 






Z\ 



^E[I,{AN)Iq{hN)] 



<l^-{fi,N T Jj.n) L'^{[o,i\i) ^iid SO the covariance structure of Z^'^'^' implies 

N-l 
(9!)^ ZZ (/i,Af'/j,w)i2([o,i],) 



N-l 

■E 

i,i=0 



^ - j + 1 

N 



2H 



N 



2H 



N 



2H' 



(4.1) 



Secondly, the square norm of the contraction fi^N®kfj,N has been computed before 
(actually its expression is obtained in the lines from formula (3.7) to formula (3.8)). 
By a simple polarization, we obtain 



7V-1 



^ 2^ {fi,N ®k fj,N, fj.N ®fe /i,Ar)L2([0,l]29-2'=) 

N-1 „l „l „l „l 

= d{H,qfa{H')^m^"-^N^^"'-^^^^J2 / / dydzdy'dz' 

ij^Q-Jo Jo Jo Jo 

(|y-z + z-j||y'-z' + z-j|)(2^'-2)'=|(y-y' + z-j||z-z' + z-j|)(2«'-2)(.-fe) 

and as in the proof of Proposition 3.1 we can find that this term has to be renor- 
malized by, ii H > j 

bH.kN^^-^"'^^^ = bH,kN^-'" 

where h^H,k = {q\f{Clf2d{H,qfa{H')'''ijl{l-x)x''"-^dx. If H e (i, |), then, 
the same quantity will be renormalized by b2.H,kN where 



b2,H.k - {qincird{H, q)MHT E i^f"'" - (^ + 1)'"" - (^ - ^)'"y 



fc=l 



while for iJ = I the renormalization is of order b3^H,kN{\ogN)'^ with bz,H,k = 
(ql)'^ {Cl)'^2d{H , q)*a(H')'^'^ (1/2). As a consequence, we find a sum whose behavior 
is well-known (it is the same as the mean square of the quadratic variations of the 
fractional Brownian motion, see e.g. [20]) and we get, for large N 



^[Tiq] 



^^iM, ifife(i,J);E[T|J 
logA^ ., „ 3 



N' 



4H-4 



X2.H, a H e (-,1) 
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where x^.h - {^tl Hha + 1 + (1/2) Y.Zi {"^k^" - (k + l)'"" - (k - lf"[ 
X2,H = (Eti Kh,i + H^2H - 1)1 {m - 3)) and a;3,H = {Y.Vi Hhj + 9/16). 

Remark 4.1. The fact that the normaUzing factor for T2q when H < 3/4 is N'^^"^ 
(in particular does not depend on H) is a teU-tale sign that normal convergence 
may occur. It is possible to show that the term in the chaos of order 2q converges, 
after its renormalization, to a Gaussian law. More precisely: a) suppose that 
H <E (^,|) and let Fm := x^ ^ VNT2q; then, as iV ^ oo, the sequence F^ 
convergence to the standard normal distribution iV(0, 1); b) suppose H = j and 

set Gat :— \ / lo at ^3 h -^29! then, as A^ — > 00, the sequence Gat convergence to 
the standard normal distribution A^(0, 1). We refer to the extended version of our 
paper on arxiv for the basic ideas of the proof. 



It is possible to give the limits of the terms T2q-2 to T2 appearing in the decom- 
position of the statistics Vat . AH these renormalizcd terms will converge to Hcrmite 
random variables of the same order as their indices (we have already proved this 
property in detail for T2 in the previous section) . This is a kind of "reproduction" 
of the Hcrmite processes through their variations. 
Theorem 4.2. 

• For every iJ G (^,1) o-nd for every k = 1, . . . ,q — 2 we have 

TV— J- 00 

where ^(2«-2fe,(2g-2fe)(// -i)+i) denotes a Hermite random variable with 

self- similarity parameter {2q — 2k){H' — 1) + 1 and 

Zk,H = d{H, qfa{H'f {{H' ~ \)k + 1)-^ (2(F' - 1) + 1)"^ 

• Moreover, if H G (|,1) then 

\im N^-^"x-]i^T2q = Z^^'^'^"-^\ inL\n). (4.3) 

Proof: RccaU that we have T2k = N'^"~^hq-2k f Z^ilo^ fi,N ®k fi,Nj , for 
k ~ 1 to 2q — 2, with /i^^r (g)fe fi^jsi given by (3.1). The first step of the proof is to 
observe that the limit of N^'^"' ~'^^'^'i~^^T2q-2k is given by 

^^ ^^ hq-2k{j2q-2k) 

with f2q-~2kiVi^ • • • , yq--k,zi, ..., Zq-k) = d{H , qfa{H'f ^^^0^ l[o__^] (2/i)l[o, j^] i^i) 
Jj^ Jj^ diK{u, 2/i) . . . diK{u, yq^k)diK{v, zi) . . . diK{v, Zq^k)\u - v\^'^"'-'^'^^dvdu. 
The argument leading to the above fact is the same as in Theorem 3.2: a the 
remainder term r2q-2k, which is defined by T2q-2k niinus ^2(/2^_2fe) converges to 
zero similarly to the term r2 in the proof of Theorem 3.2 because of the appear- 
ance of the indicator functions I7. (j/j) or li.{zi) in each of the terms that form this 
remainder. The second step of the proof is to replace diK{u,yi) by diK{j^,yi) 
and diK{v,Zi) by diK{jj,Zi) on the interval li inside the integrals du and dv. 
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This can be argued, as in the proof of Theorem 3.2, by a dominated convergence 
theorem. Therefore the term N^^^^^ ''•''^'^■'72q_2fc will have the same limit as 



N-l 



Ni^--^')i^-^)N^^-'d{H^qfa{H'f J2 l[o,-](y.)l[o,^](-.) 



q — k 

n 



d,K{^,y,)d,K{^,z,) 



i=0 



d(H,q) a(H ) ^(2-2H')(q-k) ^2H-1 ^(2-2H')k+2 



{{H' - l)k + I) {2{H' - I) + I) 

N-l q-k 

N-l q-k 



diK{^,yj)d^K{^,z,) 



Zk,H 



Ei[o.^]^(y*'^*)n 



J=l 



^^K{-,y,)^^Ki-,z^ 



'N 



N' 



Now, the last sum is a Riemann sum that converges pointwise and in i^([0, 1]^'^ ^'^) 
to Zk.nLi where 

nl 9-fc 

Li{yi,...,yq^k,zi,...,Zq^k) == / J]^ [aii^(u, yj)5ii^(u, Zj)]du 

which is the kernel of the Hermite random variable of order 2q — 2fc with self- 
similarity parameter {2q — 2k){H' — 1) + 1. The case H e (|, 1) and k = follows 
analogously. ■ 

5. Consistent estimation of H and its asymptotics 

Theorem 3.2 can be immediately applied to the statistical estimation of H. We 
note that with Vat in (1.1) and Sn defined by 

N-l 



s 



N 



ly (z[i\-z 



i=0 



we have 



1 + Vn^ N^^Sn 



(5.1) 



(5.2) 



and E [Sn] = N-^" so that H = -^fj^^. To form an estimator of H, since 

Theorem 3.2 implies that Sn evidently concentrates around its mean, we will use 

Sn in the role of an empirical mean, instead of its true mean; in other words we 

let 

^^^ log^ 

21ogiV 

We immediately get from (5.2) that 

log(l + VW) =2iJlogiV + log5'Ar ==2fi?-i/')log7V. (5.3) 

The first observation is that Hn is strongly consistent for the Hurst parameter. 
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Proposition 5.1. We have that Hn converges to H almost surely as A^ — > oo. 
Proof: Let us prove that Vn converges to zero almost surely as A^ — > oo. We 



know that Vat converges to in L (fl) as N 



and an estimation for its variance 



is given by the formula (3.6). On the other hand this sequence is stationary because 
the increments of the Hermite process are stationary. We can therefore apply a 
standard argument to obtain the almost sure convergence for discrete stationary 
sequence under condition (3.6). This argument follows from Theorem 6.2, page 
492 in [6] and it can be used exactly as in the proof of Proposition 1 in [4] . A 
direct elementary proof using the Borel-Cantelli lemma is almost as easy. The 
almost sure convergence oi H to H is then obtained immediately via (5.3). H 



Owing to the fact that Vn is of small magnitude (it converges to zero almost 
surely by the last proposition's proof), log(l + Vn) can be confused with Vn- It 

then stands to reason that (H — H) N^^^^ logiV is asymptotically Rosenblatt- 
distributed since the same holds for Vn by Theorem 3.2. Just as in that theorem, 
more is true: as we now show, the convergence to a Rosenblatt random variable 
occurs in L^ (fl). 

Proposition 5.2. There is a standard Rosenblatt random variable R with self- 
similarity parameter 2H' — 1 such that 

21 



hm E 

N—>oo 



2^2-2H' (H-H)\ogN^ C2c\'l R 







1/2 

Proof: To simplify the notation, we denote C2C^ ^j by c in this proof. Theorem 
3.2 signifies that a standard Rosenblatt r.v. R with parameter 2H' — 1 exists such 
that 

2" 



lim E 

JV-s-oo 

From (5.3) we immediately get 



N^-^"'Vn-cR 



E 

< 2E 



2N 



2-2H 



' (h -H\ log N -cR 

2' 



E 



= 0. 



Ar2-2-ff' log (1 + Vat) -ci? 



iV2-2ffVw-ci? 



2N^-^"'^ 



\VN-\og{l + VN)Y 



\VN-\0g{l + VN)Y 



o{N 



4ff'-4 



Therefore, we only need to show that E 

Using the inequality x — log(l + x) — \x — log(l — a:;)| < x^ for a; > ^ one gets 

E|Fw-log(l + V7v)|' 

< 2E \Vn - log(l + Vn^ 1v„>_ 1 + 2E \Vn - log(l + Vn)]^ h 
I' 



^Vm<-\ 



< 2EV^ + AP{Vn < 



-4E|log(l + VAr)|n 



Viv<- 



The first term is bounded above and this is immediately dealt with using the 
following lemma. 
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Lemma 5.3. For every n > 2, there is a constant Cn such that E 



\Vn\' 



< 



Proof: In the proof of Proposition 3.1 we calculated the L^ norm of the terms 
appearing in the decomposition of Vat, where Vn is a sum of multiple integrals. 
Therefore, from Proposition 5.1 of [10] we immediately get an estimate for any 
event moment of each term. Indeed, for k — q ~ 1 

E [r|"] < 1 • 3 • 5 • • • (4n - 1) {cJh N^"'-^ = c„iV(4H'-4)„ 

and for 1 < fc < g — 1 

< 1 • 3 • 5 • • • (2(2(7 - 2fc)n) {cl^^^ iV(29-2fe)(2H'-2))« 



E [n-_2k\ 



The dominant term is again the T2 term and the result follows immediately. 



Then applying the above lemma for n = 2 and using a similar proof to prove 
that E|log(l + VW)|'ly„<_i = o(iV4«'-4) and P{Vn < -^) - 0(iV«^'-*^) = 
o{N^"'^^). This leads to Proposition 5.2 ■ 

A difficulty arises when applying the above proposition for model validation 
when checking the asymptotic distribution of the estimator H: the normalization 
constant A^2-2H logA^ depends on H. While it is not always obvious that one 
may replace this instance of H' by its estimator, in our situation, because of the 
L^ (fl) convergences, this is legitimate, as the following theorem shows. 

Theorem 5.4. Let H' = 1 + {H — l)/q. There is a standard Rosenblatt random 
variable R with self-similarity parameter 2H' — 1 such that 



lim E 



2N2-2H' (^ _ ^\ log^ _ C2cl%gR 



= 0. 



Proof: By the previous proposition, it is sufficient to prove that 



hm E 

Af-5-oo 



N 



2-2H' 



N 



2-2H' 



H - HWosN 



0. 



We decompose the probability space depending on whether H is far or not from 
its mean. For a fixed value £ > which will be chosen later, it is most convenient 
to define the event 



We have 



E 
E 



D = lH>qe/2 + 2H-l\ 



j^2-2H' _ j^2-2H' ^ lH-H)\0gN 



Id 



N 



2-2H' 



N 



2-2H' 



H - HWosN 



N' 



2-2H' 



N 



2-2H' 



H -HWosN 



A + B. 
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We study A first. Introducing the shorthand notation x — max ( 2 — 2H', 2 — 2H' 
and y = min I 2 — 2H' , 2 — 2H' ) , we may write 



N' 



2-2H' 



N 



2-2H' 



^x log N _ ^y log N 



= ey '°s N L(^^v) i°g ^ -l\ <Ny (log iV) (x - y) N'^'V 



2 {log N)N'' 



H'-H' 



^2q-^ (log N)N'= 



H-H 



Thus 



A < 2q-^ E 



2q-^ E 



li^A^" 



iJ-i7 



log^TV 



N' 



x-(4-4H')- 



IdN' 



4-4_H"' 



H-H 



log^ iV 



Now choose e e (0, 2 — 2_ff'). In this case, ii lu E D, and a; = 2 — 2_ff', we get 
X ~ (4 — AH') = —X < — £. On the other hand, for w e D and x = 2 — 2H', we get 
a; - (4 - 4ii^') = 2 - 2H' - (4 - 4i/') < -e as weU. In conclusion, on D, 

X - (4 - 4iJ') < -e, 

which implies immediately 

.2 



A < N-''2q-^E 



N 



4-4H' 



H-H 



loe'^N 



this prove that A tends to as A^ — > oo, since the last expectation above is bounded 
(converges to a constant) by Proposition 5.2. Now we may study B. We are now 
operating with lu G D"^. In other words, 

H -H >1- H - qe/2. 

Still using £ < 2 — 2H\ this implies that H > H. Consequently, it is not inefficient 
to bound A^2-2H _ j\j^2-2H ^i^ove by iV^^^^ . In the same fashion, we bound 

H-H above by H . Hence we have, using Holder's inequality with the powers 
2q andp-i + (2g)"^ = 1. 

B < HlogN E Id-N- 

< HlogN Pi/P [D^] Ei/(29) rAr(2-2H')29 

From Proposition 5.2, by Chebyshev's inequality, we have 



2-2H' 



(5.4) 



eVp 



pl/P [£)C] < 



H-H 



< Cq^HN 



-{4-4H')/p 



(5.5) 



{1-H- qe/2f^P 

for some constant Cq.H depending only on q and H. Dealing with the other term 
in the upper bound for i? is a little less obvious. We must return to the definition 
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of H. By (5.3) we have 

1 + Vn = Ar2(H-H) = j^2q{H'~H') 
^ j^2q{2-2H') j^-2q{2-2H') _ 

Therefore, 



E1/(2<?) 



Ni^-^^'>i] < Ar2-2ff'Ei/(29)[i + y^] 



< 2N^-^"' . (5.6) 



Plugging (5.5) and (5.6) back into (5.4), we get 

B < 2i/c,^H(logiV)iV-M^')(2/p-i). 

Given that 2g > 4 and p is conjugate to 2q, we have p < 4/3 so that 2/p — 1 > 
1/2 > 0, which proves that B goes to as A^ — > cx). This finishes the proof of the 
theorem. ■ 

Finally we state the extension of our results to L^ (ri)-convergence. 

Theorem 5.5. The convergence in Theorem 5.4 holds in any L^ (il). In fact, the 
L^ (fl)- convergences in all other results in this paper can be replaced by L^ (fl)- 
convergences. 

Sketch of proof. We only give the outline of the proof. Lemma 5.3 works 
because, in analogy to the Gaussian case, for random variables in a fixed Wiener 
chaos of order p, existence of second moments implies existence of all moments, 
and the relations between the various moments are given using a set of constants 
which depend only on p. This Lemma can be used immediately to prove the 
extension of Proposition 3.1 that 



E 



V, 



2v 



N 



CpiVP(4ff'-4). 



Proving a new version of Theorem 3.2 with L^ (r2)-convergence can base itself on 
the above result, and requires a careful reevaluation of the various terms involved. 
Proposition 5.2 can then be extended to L^ (ri)-convergence thanks to the new 
LP (fl) versions of Theorem 3.2 and Proposition 3.1, and Theorem 5.4 follows 
easily from this new version of Proposition 3.1. H 
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